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Substantial collective flow is observed in collisions between Lead nuclei at LHC as evidenced by 
the azimuthal correlations in the transverse momentum distributions of the produced particles. Our 
calculations indicate that the Global t'l-flow, which at RHIC peaked at negative rapidities (named 
as 3rd flow component or anti-flow) , now at LHC is going to turn toward forward rapidities (to the 
same side and direction as the projectile residue). Potentially this can provide a sensitive barometer 
to estimate the pressure and transport properties of the Quark-Gluon Plasma. Our calculations also 
take into account the initial state Center of Mass rapidity fluctuations, and demonstrate that these 
are crucial for vi simulations. In order to better study the transverse momentum flow dependence 
we suggest a new "symmetrized" vf{pt) flow component; and we also propose a new method to 
disentangle Global vi flow from the contribution generated by the random fluctuations in the initial 
state. This will enhance the possibilities of studying the collective Global vi flow both at the STAR 
Beam Energy Scan program and at LHC. 

PACS numbers: 12.38.Mh, 25.75.-q, 25.75.Nq, 51.20.+d 



The first publication from the LHC heavy ion run pre- 
sented amazingly strong elliptic flow, exceeding all mea- 
surements at lower energies [1 . This indicates strong 
equilibration and thermalization at these energies in con- 
trast to expectations of increasing transparency. Just 6 
month later ALICE has also measured the vi flow [2 . 

The overall picture indicated by the first vi results is 
very similar for RHIC and ALICE/LHC; namely that Vi 
has three physical sources O |3] : 

(i) the Global collective flow correlated with the reaction 
plane of the event EP; 

(ii) the Random fluctuation flow of all Vn varieties, where 
the corresponding symmetry axes for (e.g. for Vi and v^,) 
have no correlation with the reaction plane EP, instead 
they are observed with respect to a participant plane PP 
Event by Event (EbE) [H |5] . The participant planes are 
different for the neighbouring flow harmonics; 

(iii) at high momenta or high pseudo-rapidity, (1.5 < 
1 77 1 < 4), there are strong anti-flow peaks (in opposite 
direction with respect to classical bounce-off). These ap- 
pear only at RHIC and LHC energies [2l|3], and there this 
is the strongest source of vi. These high momenta parti- 
cles are not, and probably can not be, described by fluid 
dynamical models. It seems reasonable that they are 
generated in very early (pre-equilibrium) times of the re- 
action, and such an emission is anti-correlated with the 
projectile spectator in the reaction plane due to shad- 
owing effect of the main reaction volume [6l |7]- The 
hybrid transport AMPT model provided a qualitative 
match to this Vi flow component under special assump- 
tions (switching off 'sting- melting') [6j, what basically 



means very early (pre-equilibrium?) hadronization and 
freeze out in some parts of the reaction volume. 

This article discusses the behaviour of the first (i), 
among these flow phenomena, which is the weakest at 
RHIC and LHC energies. We will also discuss, how to 
separate the Global vi-flow, from the one produced by 
random EbE fluctuations of the initial state (ii). 

Collective flow is evidenced by the radial flow and, in 
non-central collisions, by the asymmetric azimuthal dis- 
tribution around the beam axis quantified by the func- 
tions Vi{y,pt),V2{y^Pt)^ ... in the expansion 
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where y is the rapidity and pt is the transverse momen- 
tum and (j) is the azimuth angle in the transverse plane 
with respect to impact parameter vector, b . 

The observed large V2{pt) has important consequences. 
As the previously observed constituent quark number 
scaling indicates, the collective flow must have developed 
in the Quark-Gluon Plasma (QGP) phase, and the flow 
at the partonic level becomes observable after partons co- 
alesce [8 . Theoretical calculations also indicate that to 
explain the observed flow enhanced partonic interaction 
is needed over perturbative QCD predictions [9^. Thus 
the QGP is strongly interacting. At the same time theo- 
retical estimates and observations also indicate that the 
QGP is a nearly perfect fluid, with minimal shear viscos- 
ity at the phase transition point [101 E] • 

Our modeL The energy- momentum tensor density 
for a perfect fluid is T^^ = {e-\-P)u^u^ —Pg^^ ^ where P is 
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FIG. 1. (color online) Initial energy density [GeV/fm^] dis- 
tribution in the reaction plane, [x,y] for a Pb+Pb reaction 
at 1.38 + 1.38 A-TeV collision energy and impact parame- 
ter b — 0.56macc at time 4 fm/c after the first touch of the 
colliding nuclei, this is when the hydro stage begins. The 
calculations are performed according to the effective string 
rope model [12^. This tilted initial state has a flow velocity 
distribution, qualitatively shown by the arrows. The dashed 
arrows indicate the direction of the largest pressure gradient 
at this given moment. 

the local pressure, e is the local energy density, and = 
7(1, iT) is the local flow velocity. We assume the MIT Bag 
Model Equation of State during the whole calculation. 

A fluid dynamical (FD) description of the nuclear mat- 
ter is considered here for Pb+Pb collisions at 1.38 + 1.38 
A-TeV. The matter expands until it reaches freeze-out 
(FO). The FD description does not constrain the FO: an 
external condition, for example a fixed FO temperature, 
is needed. The FD model we use p^HIS] can run well 
beyond the FO, so the location of physical FO can be 
selected afterwards as a space-time hyper- surface. 

The (3+l)-dimensional, FD model pISHT? uses the 
Particle in Cell (PIC) method adapted to ultra- 
relativist ic heavy ion collisions. The numerical dissipa- 
tion of the method was analyzed recently in [16]. In this 
method, marker particles, corresponding to fixed baryon 
charge, move in an Eulerian grid. The calculation, de- 
scribing the reaction, starts from an analytic initial state 
model [12 , based on longitudinally expanding strings 
of the color-magnetic field. The produced initial state, 
shown in Fig. [T] is tilted, and, thus, the direction of the 
largest pressure gradient is pointing in the "anti-flow" 
direction, what resulted in anti-flow peaks in simulations 
for RHIC and SPS [171 EH] • However, one should not for- 
get that this initial state also has a flow velocity distribu- 
tion, which tends to further rotate it, i.e. effectively it has 
a large initial " angular momentum" , what will change the 
direction of the strongest pressure gradient with time. 



Fig. [2] shows the energy density distribution in the 
reaction plane later. One may notice that the final state 
is strongly rotated with respect to the initial one, due to 
the large initial " angular momentum" , and the direction 
of strongest transverse expansion points to O = 75/255^. 
Thus, the upward moving matter is moving now forward 
and the downward moving matter backward, in contrast 
to what happens at RHIC and SPS energies [14 . 

The fluid cells in the presented calculations were 
(0.438/m)^ for peripheral collisions, h = 0.5 — O.Tbmax- 
While initially we had 2500- 5400 fluid cells containing 
matter, this increased over 100 000 by the end of cal- 
culation. The higher energy at LHC results in a more 
explosive expansion, which leads to an explosion shell 
with decreasing central density. 

In a simplest approach we assume a constant time 
FO hypersurface. Comparing measured multiplicity b- 
dependence at LHC with our FD multiplicity, we have 
chosen tpo = 8 fm/c after the formation of the hydro 
initial state. The transition from pre FO QGP to post 
FO ideal massless pion Jiittner gas is calculated accord- 
ing to the method described in ref. [19], satisfying the 
conservation laws. In this way for each fluid cell z, we 
obtain a flow velocity = {vt'^^vl) of the gas and its 
temperature T\ 

Using the Cooper-Frye FO formula we obtain: 
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where f^{y,Pt) is the normalized momentum distribution 
for cell i] the angle (j) is taken with respect to the reaction 
plane. Then, 
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where is the flow rapidity and (j)^ is the azimuthal 
angle of the flow velocity in the transverse plane of the 
given cell i. In eq. Q we have rewritten flow 4- velocity 

in the following way: = 7^(cosh?/o7 sinhT/g, ^^J), with 
= 1/y^1 — {vl)'^. In is a Bessel 
a fb^ dx\/ — a^e~^ is a mod- 



function, and i^i(a,6) 



ifled Bessel function of the second kind. 
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FIG. 2. (color online) The energy density [GeV/fm^] distri- 
bution in the reaction plane, [x,z] for the reaction shown in 
Fig. [l]at time 12 fm/c after the formation of the hydro initial 
state. The expected physical FO point is earlier but this post 
FO configuration illustrates the flow pattern. 

Results: pt-dependence of the flow. The calcu- 
lated V2{pt) distributions are similar to the experimental 
trends. For illustration one calculated 'U2(p^)-distribution 
is presented in Fig. [3) The full curve, calculated accord- 
ing to eq. ([3|, is slightly below the experimental data. 
This can be attributed to the integral over the whole ra- 
pidity range, while the experiment is only for \r]\ < 0.8, 
and to the initial state fluctuations, as discussed below. 

As vi is an antisymmetric function of the y- 
integrated vi{pt) value must vanish. In our calculation 
this is realized to an accuracy better than 10~^^. How- 
ever, considering this obvious asymmetry, we propose 
to constuct a symmetrized function vf reversing the pt- 
direction of backward going {y < 0) particles: 



j^ceiis j'^'^dcl) p{y,^gii{y) -pt) cos( 
YT^^' lo^'d^ P(y^ sgn(y) • pt) 



(4) 



where sgn{y) extracts the sign of rapidity. The idea stems 
from Danielewicz and Odyniecz [20]. In this way we get 
a no n- vanishing Vi{pt) function, which will be also much 
less sensitive to the initial state fluctuations. 



^fi^2D{v\T\pt)h{YvlPt/T^)cos{^^ 
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where D{v,T,pt) = e'^P'^^j^^. The vf{pt) parame- 
ter calculated in this way is shown in Fig. [4] (full line). 

The ALICE team has made a detailed symmetry analy- 
sis of the low rapidity component in the acceptance range 
of the ALICE TPC [2 . They introduced new quantities 

^even/odd^^^^^^ = [vi{r],Pt) ± Vi{-r] , pt)] /2 , (6) 
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FIG. 3. The V2 parameter calculated for ideal massless 
pion Juttner gas, versus the transverse momentum, pt for 
b = 0.7b max J at t — 8 fm/c FO time. The magnitude of V2 
is comparable to the observed V2 at 40-50 % centrality (black 
stars). See text for more explanation of different curves. 

viiv, Pt) = t^r"(^,Pt) . (7) 

If we have a Global Mirror Asymmetric (MA) vi flow 
only, as in our hydrodynamical simulations, then the 
even component vanishes: 



(8) 



A non- vanishing vf^^'^{r]^pt) can only come from the Mir- 
ror Symmetric (MS) part of random fluctuation flow. 

Furthermore, if we will integrate over (pseudo-) rapid- 
ity, then, as it was discussed above, the Global Vi{pt) 
flow, correlated with the reaction plane of the event, 
is exactly zero. Thus, non-zero Vi'^^^{pt) and vl^^{pt) 
can only come from the random fluctuation in the initial 
state. Without other assumption, we should not make 
any difference for MS and MA fluctuations, and thus 
these two components should be (approximately) equal: 
^e^en^^^) = vf'^{pt). The preliminary ALICE results [2] 
clearly comfirm the simple logical sequence. So, we can 
conlcude that ^^(p^) = vf'^ipt) observed by ALICE 
collaboration come from the second vi flow source (ii), 
and can not tell us anything about Global vi flow (i). 

However, we can gain information about the pt depen- 
dence of the Global directed flow, if we repeat the same 
analysis, i.e. separation into even and odd components, 
for the Vi{y^pt) function, introduced above, eq. In- 
deed, we obtain: 
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where ^^ijiuctXPt) ^ contribution to the '^f function 
from the initial state fluctuations, which is approximately 
equal for both components, as discussed above. Thus, in 
this case we can separate the contributions from the (i) 
and (ii) sources: 

' ^ ^ (11) 
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Results: rapidity-dependence of the flow. The 

vi (y) dependence is shown in Fig. [s] (fuh hne is the ana- 
lytic solution, eq. Q). As we can see the Vi is relatively 
large in the experimental rapidity range \y\ < 0.8, reach- 
ing a peak of 26 % at ?/ = ±0.5. The most important 
change with respect to the similar simulations for RHIC 
[TF is that the Vi now peaks in "forward" direction, i.e. 
the positive peak appears now at positive rapidity. 

Qualitatively our results agree with the simulations 
performed in a microscopic transport model, namely the 
quark gluon string model [7 , where '^1(77) in "forward" 
direction was obtained. However, the authors of [7 have 
not found the reason for such a behaviour, and have qual- 
itatively attributed it to the different viscosities in the 
region with I77I < 3 and at higher preudorapidities. 
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FIG. 4. (color online) The flow parameter calculated ac- 
cording to the eq. (pi for ideal massless pion Juttner gas, 
versus the transverse momentum, pt for h — 0.7bmax, at t = 8 
fm/c FO time. See text for explanation of different curves. 

At lower energies in the same FD model calculations 
we obtained the vi peaking in the "backward" direction 
{3rd flow component) [18], of a magnitude of 5 % 
and 2 — 3 % for 158 and 65 + 65 A- GeV energy re- 
spectively. The position of the peaks also moved from 
\y\ ^ 1.5 to \y\ ^ 0.5 with energy increasing from SPS 
to RHIC. Experimentally the 3rd- flow component was 
indeed measured at these energies [3,[T7l[2j, although 
the peak values were smaller. Especially at the RHIC 
energies [3 , where the highest values were vi ^ 0.6 % 
and 0.2 % for for 62.4 + 62.4 and 200 A-GeV energy re- 
spectively. The peaks appeared at \y\ « 1 around the 
far end of the acceptance of the central TPC. Thus, at 
RHIC the vi magnitude was about 5 times smaller than 
the FD prediction. Also, the move towards the more cen- 
tral rapidities was weaker in the experiment than in FD 
calculations. 

The reason for such a disagreement is the effect of 
initial state Center of Mass (CM) rapidity fluctuations, 
which may be decisive in the case of vi, due to the sharp 
change around y = 0. 

Initial state CM rapidity fluctuations. One has 
to take into account that the CM rapidity is not exactly 
the same for all collisions, due to random fluctuations in 
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FIG. 5. (color online) The vi parameter calculated for 
ideal massless pion Juttner gas, versus the rapidity y for 
b = 0.7bmax, at t = 8 fm/c FO time. Full curve presents 
semi analytical calculations according to eq. ([2]); the vi peak 
appears at positive rapidity, in contrast to lower energy cal- 
culations and measurements. The dash-dotted and dotted 
curves present vi calculated taking into account initial CM 
rapidity fluctuations. 

the initial state, where the numbers of participant nu- 
cleons from projectile and target may not be exactly the 
same. This leads to considerable ycM fluctuations at 
large impact parameters, where the flow asymmetry is 
the strongest, while total number of participants is the 
smallest. Although, several initial state models generate 
EbE fluctuating initial states, see for example [22 , longi- 
tudinal fluctuations are not analysed up to now, neither 
theoretically nor experimentally in detail. A high accep- 
tance experiment could provide a good estimate for the 
EbE initial state rapidity, ycM [23 , which is a conserved 
quantity, i.e. it can not be changes by the system expan- 
sion, hadronization or freeze out. 

To analyze the consequences of these fluctuations, 
we assumed a Gaussian ycM distribution, centered at 
ycM = O5 with variance Sy = 1,2. 

Results can be seen in Fig. [5] dash-dotted and dotted 
lines. As expected the initial state fluctuations strongly 
reduce vi{y) at central rapidities. The resulting vi is 
still large enough to demonstrate the "rotation effect", 
discussed above, however being of the order of 1% it can 
be easily masked by the directed flow generated by the 
random fluctuations in the initial state. 

The first preliminary results from LHC [2 show '^1(77) 
at midrapidity of the order of 0.1% or less, in antiflow di- 
rection. Taking into account the error bars, the observed 
directed flow at LHC is very little, practically compared 
with 0. This might be a result of a compensation of the 
'^1(77) in "backward" (anti-flow) direction, coming from 
the random fluctuations in the inital state [22], with 
the Global directed flow in the "forward" direction, as 
predicted by our simulations. 

It is interesting to study the effects of the initial CM 
rapidity fluctuations on other observables. Fig. |6] shows 
the elliptic flow as a function of rapidity. Fluctuations 
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FIG. 6. (color online) The V2 parameter calculated for 
ideal massless pion Juttner gas, versus the rapidity y for 
b — 0.7bmax, at t = 8 fm/c FO time. Full curve presents 
semi analytical calculations according to eq. ([2]); dash-dotted 
and dotted curves present V2 calculated taking into account 
initial CM rapidity fluctuations. 



make the V2{y) peaks wider, but the magnitude is hardly 
reduced. Thus, we predict a plateau-like shape of the 
elliptic flow distribution. 

The CM rapidity fluctuations have, in principle, no 
effect of the ^/-integrated V2{pt) and vf{pt) (therefore the 
full line in Figs. [3j[4]are marked as "analytical solution 
(+ fluctuations)"). However, in the realistic simulations, 
we should not integrate over y from — oo to +oo, but only 
over the measured rapidity range, i.e. —0.8<y<0.8. 



Such a "limited range" effect is dramatic for (pt) which 
can be reduced to less than 1%, see the dashed and dotted 
lines in Fig. [4] The V2{pt) dependence is weakly affected 
(Fig. [s] dashed line). 

Interestingly, the initial ^/CM-Auctuations lead to some 
increase of the elliptic flow, V2{pt)^ putting it in a rea- 
sonable agreement with the ALICE data [1 , see Fig. [sj 
- please note that no fine-tuning was done. At the same 
time -fluctuations strongly reduce vf{pt)^ see Fig. [i] 

To summarize, our FD simulations of the LHC 
heavy ion collisions suggest that collective directed vi{y) 
flow and newly introduced vf{pt) function can and 
should be measured [23] , although these are strongly sup- 
pressed due to initial state ^/CM-Auctuations (see Figs. [Ij 
|6|. For the first time in hydrodynamical calculations we 
see is that the vi Global fiow can change the direction 
to "forward" in contrast to what happened at lower en- 
ergies. This is a result of our tilted and moving initial 
state [12], in which the effective "angular momentum" 
from the increasing beam momentum is superseding the 
expansion driven by the pressure. We have also proposed 
a new method to distinguish contributions to vi{pt) from 
Global fiow (i) and from random fiuctuations in the ini- 
tial state (ii). The method is based on vf{pt) function, 
introduced by us in this work, and consist in analyzing 
its even and odd components according to eqs. (12). 
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